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Context Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor 

spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. 

Aims. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 /jm wavelength, to study the variation of the spectral energy distributions 

(SEDs) with galacto-centric distance. 

^^ , Methods. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric 
^,1^ • distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 /jm and 

500 yum using also MlPS/Spitzer data, to derive first estimates of the dust physical conditions. 

Results. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all 

SPIRE maps (250-500 yum). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged 
^ . flux densities. The two component model, with yS fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; 
5-H ' the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is 
^ ' found to be about 60K+10K. The temperature of the cold component drops significantly from ~ 24K in the inner 2 kpc radius to 13 K beyond 
Ci , 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for/? = 1.5, averaged over the galaxy, is higher than the solar value by a 

factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33. 

Key words. Galaxies: Individual: M 33 - Galaxies: Local Group - Galaxies: evolution - Galaxies: ISM - ISM: Clouds - Stars: Formation 

1 . Introduction disk galaxy, as opposed to the nearer Magellanic Clouds, which 

are highly disturbed irregular dwarf galaxies. 

In the local universe, most of the observable matter is contained ^33 -^ ^ ^j^^ ^^^^ ^^^^-^^ galaxies; it has been ob- 

in ; in stellar objects that shape the morphology and dynamics of ^^^^^^ extensively at radio, millimeter, far-infrared (FIR), opti- 

O , their 'parent galaxy. In view of the dominance of stellar mass, ^^j^ ^^^ ^-ray wavelengths, ensuring a readily accessible multi- 

^ . a better understanding of star formation and its consequences ^^velength database. These data trace the various phases of the 

^ ■ is mandatory. There exists a large number of high spatial res- in^rstellar medium (ISM), the hot and diffuse, the warm and 

olution studies related to individual star forming regions of the ^^^^-^^ ^^ ^^jj ^^ ^j^^ ^^j^^ ^^^^^^ ^^^ ^^^-^^ pj^^^^^^ -^ ^^^■_ 

Milky Way, as well as of low linear resolution studies of exter- ^-^^ ^^ ^j^^ ^^^jj^^ component. However, submillimeter and far- 

nal galaxies. For a comprehensive view onto the physical and -^^^^^^ ^^^^ ^^ ^^-^^ ^^^ j^ ^^^ ^.^j resolutions have been 

J3 , chemical processes driving star formation and galactic evolution niissins so far 

■ it is, however, essential to combine local conditions affecting in- t . ,■ , ,- . • . ... 

dividual star formation with properties only becoming apparent ^" *e framework of the open time key project Herschel 

on global scales extended survey (HERM33ES) , we use all three instruments 

J 1, , onboard the ESA Herschel Space Observatory (Pilbratt et ^ 

At a distance of 840 kpc ( |Freedmanet al. || 1991 | ) , M33 is the 20IO) to study the dusty and gaseous ISM in M33. One focus of 

only nearby, gas rich disk galaxy that allows a coherent survey at HERM33ESison maps of the FIR co ntinuum obse rved with PACS 

high spatial resolution. It does not suffer from any distance ambi- (P pgUtsch et alJIJOlOh and SPIRE (iGriffin et alJl2010h . covering 

gutty, as studies of the Milky Way do, and it is not as inclined as ^^e entire galaxy. A second focus lies on observing diagnostic 

the Andromeda galaxy. M33 is a regular, relatively unperturbed pj^ ^„^ submillimeter cooling lines [C 11], [O i], [N 11], and H2O, 

toward a 2' x 40' strip a long the major axis with PACS and HIFI 

Send offprint requests to: C. Kramer, e-mail: kramer@iraiti . es (| de Graauw et al. || 2010|) . 

* Her.schel is an ESA space observatory with science instruments In this first HERM33ES paper, we use continuum maps cover- 
provided by European-led Principal Investigator consortia and with im- ing the full extent of M33, at 100, 160, 250, 350, and 500 yum. 
portant participation from NASA. These data are an improvement over previous data sets of M33, 
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Fig. 1. A composite 500 /im (red) and 160 //m (blue) map of 
M33. The most extended emission is traced by the longest wave- 
length map revealing the presence of the cold dust in the out- 
skirts of the galaxy. The shorter wavelengths mostly trace the 
branched spiral structure as well as distinct warm Hii regions 
and star-forming complexes (such as the ones labeled). 



obtained with ISO and 5pitee r (iHippelein et al.l2003tlHinz et al.l 
l2004t iTabatabaei et al.ll2007h . in terms of wavelength coverage 
and angular resolution. 

The total bolometric luminosity of normal galaxies is only 
about a factor of 2 larger than the total IR continuum emission 
dHauser & Dwekl l20bli). which in turn accounts for more than ~ 



98% of the emission of the ISM (dust+gas) (e.g. lMalhotraet al.l 
|200U iDale et al.l 1200 H) . Massive star formation heats the dust 
mainly via its far-ultraviolet (FUV) photons and the absorbed 
energy is then reradiated in the IR. FIR continuum fluxes are 
therefore often used as a measure o f the interstellar radiation 
field (ISRF ) (e.g. Kramer et al.ll2005l) and the star formation rate 
(SFR) (e.g. lSchuster et alll2007l) . However, a number of authors 
have suggested that half of the FIR emission or more is due to 
dust heated by a difiiise ISRF, and not directly linked to massive 
star formation (llsrael et al.lll996l:IVerlev et al.ll2009i) . 

Another disputed topic is the evidence for a massive, cold 
dust component in galaxies. The SCUBA Local Universe Galaxy 
Survey (Dunne & Bales 2001) identified a cold dust component 
at an average temperature of 21 K. A number of studies of the 
millimeter continuum emissi on of galaxies found indications 
for even lower tempe ratures jMisiriotis et al.l 120061; IWeifi et all 
120081: IlIu et al.ll2010t) . In order to estimate the amount of dust at 
temperatures below about 20 K, and to improve our understand- 
ing of the physical conditions of the big grains, well calibrated 
observations long ward of ~ ISOyum wavelength are needed. 

2. Observations 

M33 was mapped with PACS & SPIRE in parallel mode in 
two orthogonal directions, in 6.3 hours on January 7, 2010. 
Observations were executed with slow scan speed of 20"/sec, 
covering a region of about 70' x 70'. Data were taken simultane- 
ously with the PACS green and red channel, centered on 100 and 
160 jum. SPIRE observations were taken simultaneously at 250, 
350, and 500 yum. The PACS and the SPIRE data sets were both 
reduced using the Herschel Interactive Processing Environment 
(HIPE) 2.0, with in-house reduction scripts based on the two 
standard reduction pipelines. 



2.1. PACS data 

The maps are produced with "photproject", the default map 
maker of the PACS data processing pipeline, and a two-step 
masking technique. First we generate a "naive" map, i.e. not 
properly taking into account partial pixel overlaps and geometric 
deformation of the bolometer matrix, and build a mask consid- 
ering that all pixels above a given threshold do not belong to 
the sky. Then we use this mask to run the high-pass filter (HPF) 
taking into account this map. The mask helps to preserve the 
diffuse component to some extent. With new HIPE tools becom- 
ing available, we will try improving data processing to fully re- 
cover the diffuse emission in the PACS maps (cf. Fig[T]). The fi- 
nal map is built using the filtered, deglitched frames. They have 
a pixel size of 3.2" at 100 fim and 6.4" at 160 fim. The spa- 
tial resolutions of the PACS data are 6.7" x 6.9" at 100 yum and 
10.7" X 12.1" at 160/im. The pipeline processed data were di- 
vided by 1.29 in the red band and 1.09 in the green one, as 
this correction is not yet implemented in HIPE 2.oO] The rms 
noise levels of the PACS maps are 2.6 mJy pix"^ at 100 fim and 
6.9 mJy pix"^ at 160 yum. The background of the PACS maps of 
M33 shows perpendicular stripes in each scanning direction due 
to 1// noise. 

2.2. SPIRE data 



A baseline fitting algorithm dBendo et al.ll2010l) was applied to 
every scan of the maps. Next, a "naive" mapping projection was 
applied to the data and maps with pixel size of 6", 10", and 14" 
were created for the 250, 350, and 500 //m data, respectively. 
Calibration correction factors of 1.02, 1.05, and 0.94 were ap- 
plied to the 250, 350, and 500 yum maps, as this is not yet imple- 
mented in HIPE 2.0. The spatial resolutions are 18.7" x 17.5", 
26.3"x23.4", and 38.1"x35.1" at 250, 350, and 500 yum, respec- 
tively. The calibration accuracy is 15'%q The rms noise levels of 
the SPIRE maps of M33 are 14.1, 9.2, and 8 mJy/beam, at 250, 
350, 500 /im. 

3. Results 

3.1. Maps 

Figure[T] shows a composite image of the 160jum, 250 yum, and 
500yL/m PACS and SPIRE maps. All data sets show the floc- 
culent and knotted spiral arm structure, extending slightly be- 
yond 4 kpc radial distance. The PACS 160 yum map provides the 
most detailed view, thanks to its unprecedented linear resolution 
of 50 pc, allowing to resolve individual giant molecular clouds 
(GMCs) over the entire disk of M33. A large number of dis- 
tinct sources delineates the spiral arms. The pro perties of these 
sources are studied by IVerlev et al.l (1201 Oh and iBoquien et al.l 
(l2010h . The SPIRE data show a faint, difiiise disk, extending 
out to ~ 7 kpc. Outside of 8 kpc, both maps show some weak 
emission. 

Galactic cirrus is evident only in the outermost part of the 
galaxy beyond 6 kpc radial distance, showing an average con- 
tamination of the order of 2% which can go up to 8% at the very 
faint levels at 500 //m. This is still below the 15% calibration er- 
ror, which is the dominant part of the uncertainty. We did not 
correct the M33 data for Galactic Cirrus emission. 



' PACS photometer - Prime and Parallel scan mode release note. 
V.1.2, 23 February 2010 

2 SPIRE Beam Model Release Note VO.l, SPIRE Scan-Map AOT 
and Data Products, Issue 2, 21-Oct-2009 
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Fig. 2. Ratio map of the 250 //m over 500 yum maps, at 40" res- 
olution. Ellipses den ote 2, 4, 6, 8 kpc galacto-centric distances 
with i nclination 56 ° ( Regan & V ogel 1994) and position angle 
22.5° (iPaturel et al.ll2003l) Therectangle deUneates the 2' x 40' 
wide strip along the major axis, which will be mapped with HIFI 
and PACS in spectroscopy mode. Crosses mark the H n regions 
shown in Fig. 1. 



The 5(250 jum)/5'(500jum) ratio of flux densities (Fig.|2]l 
drops from about 6 in the inner spiral arms, to ~ 4 at ~ 4 kpc 
radius, continuing to less than ~ 3 at more than 6 kpc radial dis- 
tance. This drop is also seen in the radially averaged spectral 
energy distributions (Fig. [3] Table[TJ. In addition, the inner spi- 
ral arms and a couple of prominent Hii regions (cf. Fig.[T]i, out 
to about 5 kpc radius, show an enhanced ratio of ~ 6 relative to 
the inter-arm ISM, exhibiting a ratio of typically ~ 4. This shows 
that dust is mainly heated by the young massive stars rather than 
the general interstellar radiation field in M33. This is in agree- 
ment with a multi-scale study of MIPS data dTabatabaei et al.l 
l2007h . where the 160//m emission was found to be well corre- 
lated with Ho- emission. . 



3.2. Spectral energy distributions (SEDs) 

Figure|3] shows the total flux densities of M33 and radially av- 
eraged SEDs. The SEDs at different annuli were created by 
smoothing all data to a common resolution of 40", and averag- 
ing the observed flux densities in radial zones of 2 kpc width: 
n < R < n + 2 kpc with r, = 0,2, 4, 6 kpc (cf. Fig.|2li. The 
Herschel data agree in general well with the data from the liter- 
ature. The MIPS data at 160;um agree within 20% with the cor- 
responding PACS data, for all radial zones. The 100 yum PACS 
flux density, measured in the outermost annulus, is far below the 
expected value, indicating that extended, diffuse emission is at 
present lost by the data processing. We do not use these data for 
the fits. 

Figure|3};,d shows the drop of emission by almost two orders 
of magnitude between the center and the outskirts at 8 kpc ra- 
dial distance. One striking feature of the radially averaged SEDs 
is the change of the 160/250 PACS/SPIRE flux density ratio 




A[;jm] 



A[/im] 



Fig. 3. Spectral energy distributions (SEDs) of M33 at 
wavelengths between 24 //m and 500 jum. a, b: Total inte- 
grated SED of M33, com bining data of PACS & SPIRE, 
with data of MlPS/Spitzer dTabatabaei et all l2007h. IS OCAM 



( Hippeleinetal. 2003) , and IRAS dRice et al.1 Il990l) . c, d 



Radially averaged SEDs in zones of 2 kpc width. Here, we show 
the MIPS, PACS, and SPIRE data, a, c: Drawn lines show two- 
component grey body model fit results. The 100 //m PACS flux 
density, measured in the outermost zone, was not used for the fits 
(see the text). The dust emissivity index was set to /? = 2. b, d: 
Drawn lines show two-component models for/? =1.5. a-d: All 
fits are weighted by the assumed uncertainties: 7% for MIPS 24, 
70 //m (Spitzer Observers Manual v8.0), 20% for PACS, 15% 
for SPIRE. The SEDs have not been de-projected. 



(color), which drops systematically with radial distance, from 
1.7 in the inner zone, to 0.5 in the outer zone. At the same time, 
the slope of the SPIRE data turns shallower with distance, as 
already seen in Fig.|2] 

We fit simple isothermal and two-component grey body 
models to the data. Each component is described hy Sv = 
B(v, T)Ty = B{v, T)KyMd/D^, assuming optically thin emission, 
with the flux Sy, the Planck function By, the opacity Ty, the 
dust mass Mj, the distance D, and t he dust absorption coeffi- 
cient / Cy = 0.4(v/(250 GHz))^ cm^g ' dKruegel & SiebenmorgenI 
1 1 994t iKruege J2003h . p is the dust emissivity index. The fit min- 
imizes the function ;i'^ = yX(5 vnhs - S y) / AS y.nhfi)^ using th e 
Levenberg-Marquardt algorithm dBevington & Robinson|[l992h . 
with the assumed calibration error A5v,obs- The fits are con- 
ducted at 7 wavelengths using the SPIRE, PACS, and the MIPS 
data at 70 and 24-jj.m. The 24-fim data helps in constraining 
the warm component, though its emission partly stems from 
stochastically heated small grains, not only from grains in ther- 
mal equilibrium. To maintain at least two degrees of freedom 
dBevington & Robinsonlll992h in the 2-component fit, we kept 
/3 fixed to values between I and 2. These values are typically 
found in models and observations of in terstellar dust (see litera- 
ture compiled bv lDunne & Ealesll200lh . 

The fits of isothermal models do not reproduce the data well, 
the values of the reduced ^^ are very high. To a large extent, this 
is because of the 24jjm points, which clearly require a second, 
warm dust component. Two-component grey body models result 
in a much better agreement with the data. The best fitting model 
is the two-component model with/? = 1.5. The;^'^^^ values are 
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Table 1. Results of fits of one and two emission components to 
the measured spectral energy distributions (SEDs) of the MIPS, 
PACS, SPIRE data of M33 shown in Figure|3] 





Total 


(1) 


(2) 


(3) 


(4) 


Isothermal fits 












r/[K] 


29 


25 


28 


105 


37 


p 


0.5 


1.4 


0.8 


-1.8 


-0.8 


XxA 


45 


44 


45 


3 


6 


Two-component 


fits with 


-8 = 2 








r./[K] 


16 


20 


17 


14 


12 


r„/[K] 


50 


55 


51 


49 


48 


MJ[IO«Mo] 


8.0 


I.O 


2.4 


3.6 


4.0 


MJM„ 


1000 


900 


900 


1700 


5600 


X^ A 


0.41 


0.23 


0.33 


0.42 


1.8 


M^JM, 


250 


230 


180 


160 


200 


Two-component 


fits with 


^=1.5 








r./[K] 


19 


24 


20 


16 


13 


r.,/[K] 


55 


77 


57 


52 


51 


M,/[IO«Mo] 


10 


1.2 


3.0 


4.6 


4.9 


MJM„ 


500 


3800 


480 


730 


2200 


X~ i 


0.14 


0.10 


0.12 


0.20 


1.8 


M,,JM, 


200 


190 


150 


120 


160 


Two-component 


fits with 


P=l 








r./[K] 


23 


28 


24 


19 


15 


TJ[K] 


62 


- 


67 


57 


55 


M,/[106Mo] 


12 


1.6 


3.6 


5.7 


5.8 


MJM„. 


300 


- 


440 


330 


870 


■^red 


0.22 


0.6 


0.26 


0.30 


2.4 


M,,JM, 


170 


140 


120 


100 


140 


Total gas mass 












M,,J[10'M^] 


2020 


230 


440 


560 


790 



Notes. For the two-component fits, the dust emissivity index was kept 
fixed. Tc, r„. are the temperatures of the cold and warm component. Mc 
is the total cold dust mass per annulus. MclM„ is the dust mass ratio 
of both components, x^^^ is the x^ divided by the number of observed 
parameters minus the number of fitted parameters minus 1 . The lOOyum 
flux density measured in the outermost zone was not used for the fits. 
The columns give the radial annuli: Total: < i? < 8 kpc, (I): < i? < 
2kpc, (2): 2 < R < 4kpc, (3): 4 < R < 6kpc, (4): 6 < R < 8kpc. 
The last line at the bott om of the table give s the total gas masses Mgas = 
1.36 (M(H2) + M(Hi)) dGratier et al.ll20loh . 



better than 0.2 for all annuli out to 6 kpc, and 1.8 for the outer- 
most annulus. However, these values are only slightly better or 
equal to the x^^^ values of the two other two-component models. 

We find higher temperatures of the cold component for lower 
values of yS, rendering it difficult to determine both parameters at 
the same time. This degeneracy between dust temperature and 
dust emissivity is a common problem (e.g. Kramer et al. 2003]). 
Note, however, that the total mass of the cold component, is 
rather well determined. For j6-values varying between 1 and 2, 
this mass only varies by ~ 20%. The masses in each annulus 
were determined by fitting the observed SEDs. As the fitted tem- 
peratures are slightly different, the sum over the four annuli does 
not exactly agree with the fitted total cold dust mass in Table[T] 

The cold dust component dominates the mass for the galaxy 
for all annuli. Though the warm component is needed to repro- 
duce the data shortwards of ~ 100/im, its relative mass is less 
than 0.3% for all cases. Therefore, its temperature is not well 
constrained. It is found to be about 60 K+10 K. The temperature 
of the cold component is determined to an accuracy of about 3 K, 



as estimated from a Monte Carlo analysis using the observed 
data with the estimated accuracies. It drops significantly from 
~ 24 K in the inner 2 kpc radius to 13 K beyond 6 kpc radial dis- 
tance, using the best fitting model. 

Table[T] also gives the total gas masses Mgas of the entire 
galaxy, and of the elliptical annuli. These are c alculated from 
H I and CO data presented in lGratier et al.l (l2010l) . They assume 
a constant CO-to-Hi conversion factor. But note that this Xco 
factor does not strongly affect the total gas masses, as the H i is 
dominating. The gas-to-dust mass ratio for the entire galaxy, us- 
ing the best fitting dust model with /? = 1.5, is ~ 200, about 
a fa ctor of 1.5 higher than the solar value of 137 (cf. Table 
2 in iDraine et al.l f2007h . and a f actor of 2 higher than recent 
dust models for t he Milky Way dWeingartner & Draind I200l1: 
iDraine et al .1120071 ). A fact or of about 2 is expe cted, as the metal- 
licity is about half solar dMagrini et al.ll20()9l) . The gas-to-dust 
ratio foryS =1.5 varies between 200 and 120 in the different 
annuli. Within our errors this is consistent with the shallow O/H 
abundance gradient found bv lMagrini et al.l (l2009l) . The gas-to- 
dust ratios found in M3 3 are similar t o the typi cal values foun d 
in nearby ga l axies (e.g. IDraine et al1 l2007; Be ndo et al.ll2010l) . 
iBraine et al.l (1201 Ol) combine the dust and gas data of M33 to 
study the gas-to-dust ratios in more detail and derive dust cross 
sections. 
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